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EFFECT OF INITIAL MIXTURE TEMPERATURE ON FLAME SPEED OF METHANE-AIR, PROPANE- 
AIR, AND ETHYLENE-AIR MIXTURES ! 


By Gorpon L. DUGGER 


SUMMARY 


Flame speeds based on the outer edge of the shadow cast by the 
laminar Bunsen cone were determined as functions of composi- 
tion for methane-air mixtures at intial mixture temperatures 
ranging from — 182? to 842? C and for propane-air and ethylene- 
gir mixtures at initial mirture temperatures ranging from —78° 
to 844° C. The data showed that maximum flame speed in- 
creased with temperature at an increasing rate. The percentage 
change in flame speed with change in initial temperature for the 
three fuels followed the decreasing order, methane, propane, and 
ethylene. Empirical equations were determined for maximum 
flame speed as a function of initial temperature over the tempera- 
ture range covered for each fuel. — 

For each fuel it was found that, with a fixed parallel-beam 
shadowgraph system, the ratio of fame speed based on the outer 
edge of the shadow cast by the flame cone to flame speed based on 
the inner edge of the shadow was a constant, independent of 
temperature or composition. The flame speed of propane-air 
flames was independent of tube diameter from 10 to 22 milli- 
meters or streamflow Reynolds number from 1500 to 2100. 

The observed effect of temperature on flame speed for each of 
the fuels was reasonably well predicted by either the thermal 
theory as presented by Semenov or the square-root law of Tanford 
and Pease. The importance of active radicals in flame propaga- 
tion was indicated by a simple linear relation between maximum 
flame speed and equilibrum radical concentrations for all three 

fuels. Equally good correlations resulted from using either 
hydrogen-atom concentration alone or a summation of effective 
relative concentrations of hydrogen atoms, hydroxyl radicals, and 
orygen atoms and from using either flame temperatures based on 
a sodium D-line measurement for a room-temperature mixture 
or adiabatic flame temperatures. 


INTRODUCTION 


The flame speed, or normal burning velocity, of a fuel-air 
mixture is a fundamental property governing fame propaga- 
tion, which is one of the several major processes occurring 
in combustion equipment for flight propulsion. Inasmuch 
as plots of the performance data of aircraft combustors show 
that the combustion efficiency is related to the combustor- 
inlet temperature (for example, references 1 and 2), knowl- 
edge of temperature effects on fundamental combustion 


properties of the fuel-air mixture is desirable. The results 
of some earlier investigations of the effect of initial mixture 
temperature on flame speed are summarized in reference 3. 
Previous investigators have found flame speed to be related 
(approximately) to initial temperature raised to powers 
varying from 1 to 2. In the first phase of the present 
investigation, the effect of initial temperature on flame speeds 
of propane-air mixtures was determined (reference 4). It 
was found that flame speed increased with initial temperature 
at an increasing rate and that this trend was predicted by 
either a thermal theory (reference 5) or a diffusion theory 
(references 6 and 7). 

Inasmuch as the activation energy of the oxidation process 
is an important factor in the thermal-theory equations, it 
was decided to obtain flame-speed-temperature data for a 
gaseous fuel having an activation energy appreciably different 
from that of propane. Methane was selected because it has 
an activation energy of 51 Ialocalories per gram-mole (refer- 
ence 8) compared with 38 kilocalories per gram-mole for 
propane (reference 9, p. £37). 

Tt was also thought that data for ethylene, which had 
shown anomalous behavior in other correlations (reference 10) 
when compared with other hydrocarbons, would be of 
value. The activation energy of ethylene, which is approxi- 
mately 40 kilocalories per gram-mole (reference 11), is close 
to that of propane. 

The present report contains the results of a study of the 
effect of initial mixture temperature on tbe laminar flame 
speeds of methane-air mixtures over the range from — 132? 
to 342? C and of propane-air and ethylene-air mixtures over 
the range from —73? to 344? C. The flame speeds were 
computed from measurements based on the outer edge of the 
shadow cast by the Bunsen cone; the shadowgraphs were 
obtained with a parallel-beam shadowgraph system. The 
flame-speed values thus obtained at 25° C are compared 
with values obtained by other methods. 

Methane, propane, and ethylene are compared on the 
basis of relative increase in flame speed with initial temper- 
ature. Empirical equations for flame speed as & function of 
initial temperatures are presented. A comparison is also 
made among the experimental data for each fuel and the 
relative effects of temperature on flame speed predicted by 
the thermal and diffusion theories. 


1 Supersedes NACA TN 2170, “Effect of Initial fixture Temperature on Flame Speeds and Blow-Off Limits of Propane-Air Flames", by Gordon L. Dugger, 1820, and NACA TN 2314, 
* "Effect of Initial Mixture Temperature on Flame Speed of Methane-Air, Propane-Afr, and Ethylene- Air Mixtures”, by Gordon L. Dugger, 1950. 
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SYMBOLS 
The following symbols are used in this report: 
A longitudinal cross-sectional area of cone. (sq cm) 
dg number of molecules of combustible per unit volume 


of initia] mixture "a | 

constants in empirical equation for given fuel 

Cp mean specific heat 7, to T, (eal/(g)(° K)) 

specific heat at flame temperature (cal/(g)(° K)) 

D, diffusion coefficient at flame temperature (cm?/sec) 

D, diffusion coefficient. of i® radical at initial tempera- 
ture (cm?/sec) 

relative diffusion coefficient of 7“ radical with respect 
to other radicals 

D, diffusion coefficient of i radical at mean combustion- 

zone temperature (cm?/sec) 
E activation energy (cal/g-mole) . - 


in parenthesis following exp 

h height of cone (em) — : 

K constant from specific rate equation n 

ki specific rate constant for reaction between 2 radical 
and combustible material 

l slant height or length of generating curve (cm) 

n constant exponent in empirical equation for given fuel 

moles of reactants per moles of products from stoichi- 

ometric equation 

fi mole fraction or partial pressure of 7 radical in burned 
gas 

R ' gas constant (cal/(g-mole) (° K)) 

S lateral surface area (sq cm) 


base of Napierian logarithmic system raised to power 
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T absolute temperature (? K) 

T, flame temperature (° K) 

T, | initial mixture temperature (° K) 

u flame speed (cm/sec) 

M thermal conductivity at flame temperature 
(cal/(em?) (sec) (^? K/em)) 

Pr density at flame temperature (g/cm?) 

po density mixture at initial temperature (g/cm?) 


EXPERIMENTAL PROCEDURE 


The. flow system, high-temperature burner, and optical 


system are diagrammatically illustrated in figure 1. For 
low-temperature data, the same apparatus was used, except 
that the low-temperature burner illustrated in figure 2 was 
substituted for the high-temperature burner at section X-X 
in figure 1. 
(or cooled), and burned above a vertical tube. 

Fuels.—The minimum purities claimed by the suppliers 
of the methane and ethylene were 99.0 and 99.5 percent, 
respectively. The propane had a minimum purity of 95.0 
percent, the principal impurities being ethane and isobutane. 
Laboratory service air containing approximately 0.3 percent 
water by weight was used. 


Metering system.—The fuel and air flows were metered 
by sets of critical-flow orifices (reference 12). For each ori- 
fice, when the upstream pressure was measured with a 100- 
inch mercury manometer, the ratio of highest measurable 
flow rate to lowest was 2:1; for a set of six orifices, the over- 
all ratio of measurable flow rates was 64:1. "The six orifices 


The fuel and air were metered, mixed, heated - 
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FIGURE 1.—Diagrammatioc sketch of experimental apparatus, showing high-temperature burner. X-X, Section at which high- and low-temperature burners are exeliunged. 
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were mounted in parallel in & manifold, five of them being 
capped, while the one in use was connected by a flexible hose 
to the mixer. Ruby bearing-jewels of appropriate sizes were 
used as the orifice plates. Upstream temperatures were meas- 
ured by iron-constantan thermocouples. 

High-temperature burners.— The burners used for the 
data obtained at and above room temperature were straight 
lengths of brass tubing of the following inside diameters: 
for methane, 15.7 millimeters; for propane, 10.2, 15.7, and 
229 millimeters; and for ethylene, 6.3 millimeters. Burner 
lengths of at least 70 diameters were used, which insured 
laminar flow. Each burner tube was wrapped with asbestos- 
covered resistance wire and insulated to permit control of 
the tube-wall temperature. A 76-millimeter-diameter brass 
collar was silver-soldered to the lip of each tube to give & 
fiat horizontal surface above the insulation. <A small 
resistance heater was fastened to the under side of the collar 
for the methane and ethylene studies; in the earlier propane 
studies the lip temperature was established by the heating 
due to the flame itself. The heater section, which consisted 
of a 150-centimeter length of 9.5-millimeter-inside-diameter 
stainless-steel tubing, was attached by a union to the base 
of the burner tube. The heater tube was also wrapped with 
resistance wire and insulated. 

Temperatures were measured by iron-constantan thermo- 
couples installed at the burner-tube inlet and at the burner- 
tube port for both wall and gas temperature. The thermo- 
couple for measuring the gas temperature at the port 
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A Ceramic nozzke, 12.7 or &3 mm 
LD, with thermocouple 

B Steel flange, 20 cm diameter 
by 0.64 cm thick 

C Brass nozzle (Mache-Hebra type} 

D /00-Mesh calming screen 

E Bross collar 

F Neoprene O-ring 

-7 G Sereen-support ri 

H Brass straight sachon 

| inlet tube 

J /6 Holes equally spaced 
circumferentially 

K Manifold 

L Air jacket (sfee/ can) 

M Dewar vessel 


(thermocouple C) was of the aspirating type (reference 13) 
and was placed directly over the center of the port between 
runs. The gas temperatures at the port and at the burner 
inlet (thermocouple O) were maintained within 10? C of 
each other. 

Low-temperature burner.—For the low-temperature data, 
a nozzle-type burner was used in order to obtain laminar 
flow in a short over-all burner length, so that the whole 
burner assembly could be submerged in an 18-centimeter- 
diameter by 33-centimeter-deep Dewar vessel. The burner 
assembly is illustrated in figure 2. 

The burner proper was made of brass and had an over-all 
length of 30.5 centimeters. The fuel-air mixture was fed 
through the inlet tube | to the manifold K, from which it 
entered the burner by means of 16 small holes J, equally 
spaced circumferentially about the base. It then passed 
through the straight section H and through four calming 
screens D, which were supported by brass rmgs G and 
sealed by neoprene O-rings F. The mixture was then 
accelerated by a Mache-Hebra type brass nozzle C (ref- 
erence 14) and further accelerated by a small ceramic 
nozzle A, which was cemented to the brass nozzle. 

The ceramic nozzle was added to the original brass nozzle 
when it was found that such an arrangement aided in stabi- 
lizing the fame. The improved stability was presumably 
due to the reduced thermal gradient between flame and 
nozzle, since stable, regularly shaped flames were obtained 
only after the flame had been allowed to heat the ceramic 
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FIGURE 2.—Detad of low-temperature burner. 
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nozzle to & temperature 60? to 100? C higher than the gas 
temperature. Tests made with a propane-air mixture at 
room temperature indicated that this difference between 
the temperature of the short ceramic nozzle and the gas 
temperature did not appreciably affect flame speed. For 
methane and propane, the ceramic nozzle had a throat 
diameter of 12.7 millimeters; for ethylene, 6.3 millimeters. 
These nozzles were turned from round stock by cutting an 
axial hole to the desired throat diameter and then rounding 
the upstream edge of the hole to.approximately a 1.6- 
millimeter radius. The nozzle disks were 3.2 millimeters 
thick. 

In figure 3, a velocity profile is shown for the 12.7-milli- 
meter ceramic nozzle (on the brass nozzle) for an average 
air-flow velocity (volumetric flow rate divided by nozzle- 


throat area) of 124 centimeters per second at 25° C. The 


profile was obtained by means of s hot-wire-anemometer 
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FIGURE 3.— Velocity proflla above 12.7-millimeter ceramic nozzle cemented to brass nozzle, 
Air temperature, 25° C; Reynolds number, 1020. 


probe, whieh was ealibrated the same day by traverses of 
fully developed laminar flow (reference 15) in a 25-millimeter- 
inside-diameter tube. The calibration curve was established 
by assuming the local velocity at the axis of the 1-inch tube 
to be twice the average velocity and the average velacity to 
be equal to the local velocity at a distance from the tube 
axis equal to 0.707 times the tube radius. This condition 
was experimentally confirmed within +5 percent, which is 
believed to be the accuracy of the measurements. The 
average velocity computed from the anemometer readings 
was 131 centimeters per second, assuming a linear velocity 
gradient from the outermost data points to the nozzle wall. 
While the central portion of the velocity profile is fairly flat, 
the local velocities over this portion are of the order of 20 
percent higher than the average velocity (volumetric flow 
rate divided by nozzle-throat area). Since flame-speed val- 
ues based on a local velocity, which was assumed to be equal 
to the average velocity, would thus have been in great error, 
it was decided to use the total-area method of flame-speed 
measurement (described heremafter) for nozzle flames as 


well as tube flames. While the two-stage nozzle referred to 
herein is an unusual case, it is believed that an appreciable 
error would result from ignoring the boundary-layer effect 
for any small-diameter nozzle. 

For the data obtained at —73° C, the burner as deseribed 
was submerged in a bath of dry ice in acetone and was 
supported on the lip of the Dewar vessel M by the steel 
flange B, which was screwed onto the brass nozzle. For the 
data at — 132° C, for which liquid nitrogen was used as the 
coolant, it was found necessary to jacket ihe lower thrce- 
quarters of the burner with air to prevent overcooling and 
consequent condensation of oxygen from the primary air. 
This jacketing was accomplished by bolting a steel can L to 
a flange soldered to the brass nozzle. The desired nozzle- 
outlet temperature was then obtaimed by varying the level 
of liquid nitrogen with respect to the top of the air jacket. 

The gas temperature was measured between runs by a 
bare 28-gage copper-constantan thermocouple, which had a 
soft-solder head of approximately 1.5-millimeter diameter 
at the junction. It was supported in a piece of pyrex tubing 
with a right-angle bend such that the thermocouple leads 
extended downward (along the axis of the gas stream) for a 
distance of 6 centimeters. The wall temperature of the ce- 
ramic nozzle was measured by another copper-eonstantarr 
thermocouple. For the —73? C data, it was believed that 
errors in the gas-temperature readings duc to radiation and 
lead-conduction losses were of the order of +1 percent. For 
the —132? C data, errors due to radiation and lead conduc- 
tion when the thermocouple was placed in the nozzle throat 
were greater. A correction was estimated by the following 
method:. Gas-temperature readings were taken with the 
junction at depths of 0.3 and 5.0 centimeters below the top 
surface of the ceramic nozzle. Of these two temperature 
readings, the 5-centimeter reading was subject to smaller 
lead-conduction and radiation errors because the junction 
was well inside the cold nozzle. With the two gas-tempera- 
ture readings and the ceramic-nozzle temperature known for 
flow without burning, an estimate was made of the correc- 
tion to be applied to the 5-centimeter reading in order to 
obtain the correct temperature of the gas leaving the nozzle 
throat (initial mixture temperature). It was assumed that 
the lead-conduction error was negligible; the radiation and 
wall-to-gas heat-transfer corrections (estimated from refer- 
ence 16) totaled 5? C. ‘The initial mixture temperature was 
therefore taken to be the 5-centimeter reading plus 5° C. 

Determination of flame speed.—Shadowgraphs of the 
flames were made by a parallel-beam system as shown 
schematically in figure 1. Flame speeds were determined 
from the shadowgraphs by the. total-area method, wherein 
the average normal flame speed is equal to the volume rate 
of flow of the unburned mixture divided by the surface area 
of the cone formed by the combustion zone. This surface 
area, was determined by assuming that the flame surface 
can be approximated by the relation for conical surfaces 
of revolution; thus, 


S=qaA l/h = (1) 
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A typical shadowgraph is shown in figure 4. The area A 
was determined by measuring a magnified image of the 
shadow cone with a planimeter; the magnified image was 
obtained by tracing the outer edge of the shadow from the 
projected image of the shadowgraph. Because of inaccu- 
racies in tracing and in using a planimeter, & scatter in the 
computed points resulted, especially for the tall slender 
ethylene flames. The precision eat the higher temperatures 
was 4 percent for ethylene. This scatter was diminished 
by: (1) determining the average value of A/h for all the 
computed points at a given temperature and flow rate; 
(2) recomputing flame-speed values for one of the points on 
the basis of the average .4/h; and (3) computing the rest of 
the flame-speed values from the already corrected point by 
a ratio of cone heights according to the equation t2—=m,A;/he. 
This procedure essentially amounts to the method used in 
reference 17. 


RESULTS AND DISCUSSION 


Flame speeds based on outer edge of shadow.—.As indi- 
cated in references 11 and 18, the maximum density gradi- 
ent is more nearly represented by the outer edge of the 
shadow than by the inner edge because of the manner in 
which the shadow of a flame cone is produced. If the locus 
of maximum density gradient is taken to correspond to the 
flame front, absolute values of flame speed should be based 
on the outer edge of the shadow. In order to obtain the 
relation between flame-speed values based on the inner edge 
of the shadow (which is more sharply defined) and values 
based on the outer edge of the shadow, a group of points 
was selected at each of six temperatures from room temper- 
ature to 344° C for each fuel and was computed on both 
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bases. Each group consisted of one point representing the 
leanest mixture studied, one point near the composition for 
maximum flame speed, and one for the richest mixture 
studied. Plots of 18 points so selected for each fuel revealed 
no trends for the relating factor with respect to either mix- 
ture composition or temperature; the arithmetic mean value 
of the factors for all 18 points for & given fuel was therefore 
considered to be the best value over the ranges of composi- 
tion and temperature studied. Inasmuch as the same con- 
stant factor applies for all the points for & given fuel, the 
relative effects of temperature on flame speed would be the 
same for flame speeds calculated from either the inner or the 
outer edge of the shadow. In the early part of the present 
investigation, values based on the inside edge of the shadow 
were used because the required measurements were more 
easily and precisely made on this basis. 

The arithmetic mean values of the factors required to 
convert inner-edge values to outer-edge values of flame 
speed for the three fuels are presented in table I. The 
meximum flame speed at 25? C was obtained for each fuel 
by interpolation from figure 6 in order to compare the value 
with values obtained by other investigators who used shadow- 
cone, visible-cone, and tube measurements (references 11 
and 19 to 24). "The agreement among outer-edge shadow- 
cone, visible-cone, and tube measurements for methane and 
propane is good. The poorer agreement among various 
methods and investigators for ethylene is probably due, in 
part, to the greater uncertainties in measurements of the 
tall slender fíame shapes encountered with ethylene in both 
Bunsen burner and tube methods. 


TABLE I—FACTORS FOR CONVERTING FLAME SPEEDS 
BASED ON INNER EDGE OF SHADOW TO OUTER EDGE 
AND COMPARISON OF 25° C VALUES WITH VISIBLE 
BUNSEN CONE AND TUBE MEASUREMENTS 





Conversion factor Maximum nune at 25° E 
Fuel 
MIethane...........- : 
45. 4 
Propane.............- {5.5 i .2 í 
16 
c 76.4 
Ethylene__________.. . 849 027 l t 11.6 64. 0 f 60.0 08. 3 
i 72.0 | 


« Obtained from fig. 5 by Interpolation. 

+ Based on total area of flame traveling at uniform velocity inside a pyrex tube (reference 18) 

« Related to outer-edge value by conversion factor. 

€ Total-area method, inner boundary of luminous zone (reference m 

* Total-area method, center of Jaminous rone (present investigation}. 

f Total-area value ex lated to diametral axis of flanie; value of 45.0 cm/sec obtained 
by stroboscopieally HI ated particle method (reference 20). 

z Based on ares of upper part of shadow cone only (reference 21). 

L Based on area of cone frustrum (reference 22). 

$ Total-area method (reference 23). 

i Mean value of a number of measurements by angie method at different cone positions 
(reference L1). 


Ín figure 5, flame speed is plotted against equivalence ratio 
(fraction of stoichiometric fuel-air ratio) at several initial 
mixture temperatures for each fuel. It ts believed that no 
preflame reaction has occurred during the short preheating 
period (less than 2 sec) at these temperatures (reference 9, 
pp. 406, 411). The stream-flow Reynolds number was ap- 
proximately 1500 for methane-air mixtures and 2000 for 


110 ` REPORT 10601—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Flame speed, cm/sec 
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(8) Methane-alr flames. 
(b} Propane-air flames. 


FievRE 5.—Flame speed as function of equivalence ratio (fraction of stoichiometric fuel-air 
ratio} at various initial mixture temperatures. 


ethylene-air mixtures. The propane-air data were taken at 
two Reynolds numbers, 1500 and 2100, and are plotted 
together, showing that Reynolds number has no appreciable 
effect on flame speed in this range. Preliminary propane- 
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TABLE IL—COMPARISON OF PROPANE-AIR FLAME-SPEED 
VALUES OBTAINED WITH VARIOUS TUBE DIAMETERS 


Initíal EOD eas Tube diameter 
ture (°C) (mm) 


25 


10. 2 
18.7 
2212 
10.2 
15.7 
22.2 
3 
7 
2 


BP 





air data obtained by photographing the luminous flame cone 
above 10.2-, 15.7-, and 22.2-millimeter tubes showed the tubo 
diameter to have little effect on maximum flame speed. A 
few comparative values are given in table II. 

Each of the curves in figure 5 shows that the maximum . 
flame speed for a given temperature occurs at a mixture 
composition richer than stoichiometric (equivalence ratio 
greater than 1.0). These flame-speed maximums from fig- 
ure 5 are plotted against initial mixture temperature in 
figure 6 to show the effect of temperature on maximum flame 
speed for the three fuels. It appears that the curves in figure 
6 could be extrapolated to zero flame speed at zero degrees 
absolute. 
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FIGURE 6.—Effect of initi] mixture temperature on maximum flame speed. 
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It was found that the effect of initia] mixture temperature 
(T, °K} on maximum flame speed (u, cm/sec) could be 
represented by empirical equations of the type 


u=b+c¢T,* 


where 5, c, and n are constants for a given fuel. The equa- 
tions, which were determined by picking an integer for b 
which gave a straight line on a logarithmic plot of (u— 5) 
against 7, and then determining c and n by the method of 
least squares, are: 

For methane at 7,—141? to 615? K, 


u=8-+0.000160 T$ 
for propane at 75,—200? to 616? K, 

u=10+0.000342 Tp” 
and for ethylene at 7,—200? to 617° K, 

u=10+0.00259 Ta 7 


The logarithmic plots of the equations together with the 
data points are presented in figure 7 

In figure 8, the percentage increase in maximum flame 
speed, based on the interpolated value of flame speed at 
25° C, is plotted against temperature. The fiame speeds 
of these three fuels are again seen to be affected by tempera- 
ture in the decreasing order: methane, propane, and 
ethylene. 

Comparison of experimental data with relative values 
predicted by theoretical equations.— Two proposed mecha- 
isms for the propagation of flame are considered: One is 
based primarily on the conduction of heat from the com- 
bustion zone into the unburned gases and the other is based 
on the diffusion of active radicals from the combustion zone 
into the unburned gases. Both these theories are used to 
predict the relative effect of temperature on flame speed. 

It has been pointed out that early heat theories (for 
example, reference 9, p. 113, and reference 5), which assumed 
that (s) reaction begins at the self-ignition temperature, 
and that (b) the reaction rate is constant between the self- 


ignition temperature and the flame temperature, are inade- 
quate. The concept of ignition temperature has no meaning 
apart from an autoignition experiment, which requires a 
certain induction period. Actually, reaction rate contin- 
uously increases with temperature because it is a function 
of exp( — E/R T). 

Semenov (reference 5) derived an approximate equation 
that involves only the initia] mixture temperature, the flame 
temperature, and the physical properties at these tempera- 
tures. If the controlling step is a monomolecular reaction, 
the appropriate final equation is 

ary 


2X, Kc, , T )G ) 
pes c,p D 


and for a bimolecular reaction, the final equation is 
hee? (a) (Sav), a) Ce) 
pote? c,pD (D —Iy 
(3) 


Equations (2) and (3) take into account the cases where 
(a) specific heat and thermal conductivity vary with tem- 
perature, (b) the number of molecules changes during 
reaction, and (c) the diffusion: coefficient D, does not equal 
the coefficient of temperature conductivity (A/epp); Case (c) 
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FIGtRE 7.—Emplirical equations for maximum flame speed as function o£ Initial mixture 
temperature. 
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means that although the sum of the thermal and chemical 
energies will be the same for the over-all process, the sum of 
the energies wil vary within the combustion zone. For 
example, if D,(X/c,p), the thermal energy supplied by 
heat conduction upstream of the flame zone will be less than 
the chemical energy conducted away by diffusion. Thus, 
diffusion enters only as it affects the energy balance. 

By elimination from equations (2) &nd (3) of the terms 
independent of temperature, by the substitution, as approxi- 
mate relations for the temperature-dependent terms, of those 
relations determined for air, and by combining terms, the 
monomolecular equation can be reduced to the form 


uoc AJ T, 27" 3.88 p ee E" 


and the bimolecular equation may be reduced to the form 


uc [ret ss SEAT) (5) 


(physical properties for air were estimated from reference 16, 
pp. 391-411 and related to temperature as follows: Aoc 79; 
Eps 19-9 6 oc TO": Docufpoc TE, where u is viscosity: 
pxl'; and aycT 1) These reduced forms of the equa- 
tions ane be used to estimate the relative effect of tempera- 
ture on flame speed (thermal theory) for hydrocarbons burn- 
ing with air, provided that the relations between the physical 
properties and the temperature for the mixture are reasonably 
near those for air. 


A diffusion mechanism of flame propagation may also be : 
considered for the prediction of the effect of temperature on ` 


flame speed. Tanford and Pease (references 6 and 7) have 
proposed that the rate of diffusion of hydrogen atoms up- 


stream of the flame front determines the rate of flame prop- | 


agation. This concept, generalized to include other radicals, 
which might be important for other systems, is expressed by 


f 
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Change in maximum flame speed, percent 
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FIGURE 8.—Percentage change of maximum flame speed (based on 25° ' C) with initial mixture 
perature, 
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what lias been called the square-root law of burning velocity 


(reference 7): | 
uS 25 PDk: iL 55 (6) 


As a first approximation, it may be assumed that only 7), 
L, and p, are appreciably temperature-dependent. In ob- 
taining the temperature dependence of D,, it is necessary to 
revert to the derivation of equation (5) (reference 6) where 
D, was substituted for Dm, the diffusion coefficient at the 
mean combustion-zone temperature 7, by the relation 


EFFECT OF INITIAL MIXTURE TEMPERATURE OF FLAME SPEED OF HYDROCARBON-AIR MIXTURES 


iv pn Fe 
Dera DTP 


where 7, was assumed to be 0.7 J; Now if D, is assumed 
proportional to 74! ^ (the exponent which was used in equa- 
tions (4) and (5)), 1t follows that 


dv € Tg 


Dix C 09.7 (0.7 Ty (0.7 t Tp) Vo T 





2 
Since Loc Tp! I and DrD X Dom where D. 
f 


is the relative diffusion coefficient of the radical with respect 
to the other radicals, the temperature dependence of flame 
speed predicted by equation (6) is 


eX (SP Dar) Ts " 


Comparisons of the relative effects of initial temperature 
on maximum flame speed as predieted by the reduced equa- 
tions (4), (5), and (7) for methane, propane, and ethylene are 
presented in figure 9. The values used in plotting the curves in 
figure 9 (and fig. 10, discussed later) are presented in table III. 
The flame temperatures used in figure 9 were based on 
sodium D-line measurements of flame temperatures of mix- 
tures at room temperature (reference 25). 'Phe changes in 
flame temperature with initial temperature were assumed to 
be the same as for the theoretical adiabatic flame temperature 
(reference 26); that is, the difference between computed adia- 
batic flame temperature and sodium D-lme temperature was 
taken to be constant for a given mixture. The equilibrium 
radical concentrations were computed by the graphical 
method of reference 27. 
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From figure 9 and table III it may be seen that both the 
thermal-theory equations and the square-root law predicted 
(relative to the 25? C values) flame speeds within approxi- 
mately 20 percent for the initial temperatures and gases 
studied. The comparative agreement between the different 
theoretical curves and the experimental curve varies for the 
three fuels. In general, for elevated temperatures the ther- 
mal-theory curves are high and the square-root-law curves 
are low; whereas, for low temperatures, the values predicted 
by the two theories are substantially the same. 

It must be emphasized that even qualitative observations 
regarding the ability of one or the other of the theories to 
predict the effect of initial temperature on flame speed are 
subject to the assumptions made in deriving both the original 
equations and the reduced forms and to the choices of values 
for flame temperatures and activation energies. The effects 
of using different values of Z,and E may be illustrated by the 
case of propane at an initial mixture temperature of 343° C. 
In table III the flame-speed values predicted by the thermal- 
monomolecular, thermal-bimolecular, and square-root-law 
equations are listed as 139.0, 150.8, and 126.9 centimeters 
per second, respectively; whereas, 

(1) Using the adiabatie flame temperatures rather than 
the temperatures based on the sodium D-line measurement 
resulted in predicted values of 137.7, 149.1, and 148.0 
centimeters per second, respectively. 

(2) Using #=25,000 instead of 38,000 calories per gram- 
mole resulted in predicted values of 125.6, 136.3, and 126.9 
centimeters per second, respectively. 

(3) Combining items (1) and (2) by using adiabatic fame 
temperatures and E=25,000 calories per gram-mole resulted 
in predicted values of 124.6, 135.3, and 148 centimeters per 
second, respectively. 


TABLE IIL—PREDICTED RELATIVE FLAME SPEEDS 








































Flame temperature | Equilibrium radical pressure (total pressure of 1 atm) A RUE IDEE: mid 
Initial e (stm) speed at 25° C (cm/sec) pnr 
Equivalence ratio s kai s 
lm ee Da (edia- | pm (sodium Bor (sodium| po (sodium 
bate c | D-line 4 batic) D-line} D-line) PO -Hfne) 
[| 
L| 
2053 2203 e. Z 0x10 0. 293 X101 0. 2892x107 0. 003x104 10. 5 10.1 11.9 13. 4 
209) 2045 .310 . 298 . 107 005 17. 5 17.1 18.8 19. 2 
2153 $ 2108 . 520 - 905 . 635 E 34.0 ai Q 34. 0 734.0 
2159 2114 . 908 . 108 . 065 23.9 36. 0 35.5 35.5 
pee 1.22 (10.3 vol. percent) : 2197 2152 .6:0 .517 - 890 š 30 50. 9 41.3 48.0 
: 2231 2186 810 . 630 L 13 66. 0 08.2 60.0 62.6 
2206 2221 . 930 . 770 1.15 85.6 89.7 71.6 80. 0 
2203 2253 1.16 .9015 1.80 084 IOT. T 114. 3 90. 5 95. 0 
2345 2300 L 48 1.17 2. 14 * 14S. 2 160. 9 118.3 133.0 
T3 2204 2119 .100 - 299 1.17 " 23.8 22.2 23.0 23. 5 
2253 2193 - 599 44) 1, 62 " 10. 2 10. 2 10.2 F 10.2 
| 3| 2| 2m S |2 jm |: &i Ei i| &i 
-= & + 
Propane...... 1.49 (4.2 vol. percent) 2417 2262 | 860 Le 2.38 ] 10.2 72.5 67.8 60.2 
20 2345 2290) 1.06 . 197 2. 80 á 86.8 91.0 82. 6 81.0 
2211 2316 1.16 . 855 3. 26 í 105.7 112.3 ge. 105. 0 
2412 2351 1. 43 LOZ £.10 à 138.0 150.8 126. 9 141.2 
2302 A108 IL. 13 . 629 -9t 33.8 35.0 37. 5 35. & 
2352 2248 E £8 . 855 L 28 061.0 61.0 64.0 :6£0 
2356 2252 LS 375 1.32 O64 66. 9 67.0 66. 8 66. 6 
Ethylene 2256 2252 1.72 1.02 L 64 $8.3 $9.8 86.0 86.0 
Peery asta es 2413 2209 | L% L I7 1. Oa : HL3 114.8 106. L 105.7 
24-40 2336 2.23 La 2.33 š 137.6 144.0 128.0 126.0 
f 2185 2361 2.20 L53 2. 70 í 164 1 T:6. 1 151.0 153. 0 
2504 2100 3. 00 L83 3.39 è 218. 2 236.8 189. 8 195.0 
cm EE Mc NEM MENU L EM al EMLL LO Lii Aul 
s Fraction. of stoichiometric fuel-air ratio at which maximum flame speed occurred. r Computed from equation (7). 
+ Computed by method of ds bud ty k Reference 8. 
« Computed by method of reference i Computed using same difference between adiabatic and sodium D-line values as was 
4 Computed by subtracting from witabatte flame temperature a constant equa! to difference found for 10.0-percent methane. 
between adiabatic and sodium D-ine values (reference 24) at 25? C. i Obtained from fig. 6 by Interpolation. 


« Computed from equation (4). 
f Computed from equation CR 


k Reference 9, p. 137. 
1 Reference I 11. 
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The difference between adiabatic and sodium D-line flame 
temperatures can be seen to increase from methane to pro- 
pane to ethylene. It might be argued that since this differ- 
ence increased from fuel to fuel in the same order as the 
sodium D-line flame temperatures increased at the reference 
initial temperature of 25° C, the difference should not be 
constant for a given fuel, but should increase as the flame 
temperature increases. In other words, the sodium D-line 
temperature of 2400° K for ethylene at an initial temperature 
of 344? C might be too high. On this premise, the predic- 
tions of equations (4), (5), and (7) were recomputed for 
ethylene for a flame temperature of 2350° K at the initial 
temperature of 344° C and were found to give predicted 
flame-speed values of 196.6, 214, and 170.6 centimeters per 
second, respectively. These predicted values are well within 
20 percent of the experimental value of 196.0 and compare 
with the predicted values of 218.2, 236.8, and 189.8 centi- 
meters per second listed in table III. 

It is assumed herein that the variation of B,, & term that 
appears in the square-root law to allow for radical recom- 
bination, with initial temperature is relatively small and 
would not appreciably affect the curves in figure 6. The 
temperature dependence of k, the rate constant which 
appears in the square-root law for interaction between the 
2" radical and the combustible material, has been neglected. 
Constant values of k; were used for hydrogen, carbon mon- 
oxide, and methane flames in reference 7 where the variations 
in flame temperature are of the same order of magnitude as 
those considered herein. If the temperature dependence of 
‘ky is represented by exp(—#/RT;,), the temperature effect 
on flame speed due to E, will be small, provided that the 
activation energy for the radical-hydrocarbon reaction is 
small. For example, with an activation energy of 7 kilo- 
calories per gram-mole (reference 7}, the relative flame 
speeds predicted by the square-root law for methane, propane, 
and ethylene would be raised approximately 8, 6, and 5 
percent, respectively, at an initial temperature of 344° C 
and less at lower initial temperatures. 

The relative merits of the theories could be better assessed 
if reliable data for the effect of pressure on flame speed were 
available. The thermal-theory equations presented in ref- 
erence 5 predict that for a monomolecular reaction, the flame 
speed should be inversely proportional to the square root of 
the absolute pressure; whereas, for a bimolecular reaction, 
the flame speed should be independent of pressure. The 
square-root law predicts that flame speed should vary approxi- 


mately inversely as thefourth root of the pressure(reference 7). 


If one of the theories could be shown to be the more 
consistent in predicting both pressure and temperature 
effects, it would be strengthened considerably. Although 
some pressure data exist with which this precept might be 
tested, conflicting trends for the variation of flame speed 
with pressure for a given fuel have been reported by different 
investigators, or by the same investigator using different 
experimental methods, which indicates that some of these 
data were influenced by the experimental apparatus. For 
example, in reference 28 the flame speed of ethylene measured 
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Fiaurs 10.—Correlations between maximum flame speed and calculated equilibrium radical 
concentration. 


by & Bunsen burner method was found to be inversely propor- 
tional to the fourth root of pressure; whereas, according to 
reference 29, the flame speed of ethylene measured by a soap- 
bubble method was found by the same investigator to be 
independent of pressure. 


EFFECT OF INITIAL MIXTURE TEMPERATURE OF FLAME SPEED OF HYDROCARBON-AIR MIXTURES 


Correlation of fame speeds with active-radical concen- 
trations.—A correlation between flame speed and ?75,D,, 
F 


for ethylene-air mixtures of various compositions at room tem- 
perature is found in reference 11. In figure 10 (a), the 
maximum flame speeds from figure 4 are plotted against 
$1pD,. computed from radical pressures at temperatures 


based on sodium D-line measurements in table I; in figure 
10 (b), maximum flame speeds are plotted against computed 
hydrogen-atom concentrations at these same temperatures; 
and in figure 10 (c) maximum flame speeds are plotted against 
hydrogen-atom concentrations at theoretical adiabatic flame 
temperatures. In all cases, straight-line correlations are 
obtained. The equations of these straight lines are given 
on the figures. 

The curves of figure 10 (c) indicate that, at least for the 
three fuels studied, the maximum flame speed of & gaseous 
fuel may be accurately determined over a considerable range 
of initial mixture temperature when the experimentel values 
at two temperatures are known and sufficient data exist to 
compute the adiabatic flame temperatures and, hence, the 
equilibrium hydrogen-atom concentrations. It is recognized 
that the slopes of the correlation lmes in figure 10 are much 
steeper than those shown in the correlations in reference 10. 
In considering this difference, it must be remembered that, 
for a given correlation line in figure 10, initial mixture tem- 
perature varied while equivalence ratio (composition) and 
fuel type were fixed; whereas in reference 10 fuel type varied 
while initial temperature was fixed and equivalence ratios 
were essentially the same. The data reported herein indicate 
that flame speeds for a given fuel vary from fivefold to tenfold 
over the range of temperatures studied, whereas in reference 
10 a range of only twofold was covered in the flame speeds 
for different fuels at the same initial temperature. However, 
the flame-temperature and the radical-concentration ranges 
are of the same order in both studies. The reasons for the 
difference in slope may, in part, be attributed to the temper- 
ature dependency of the other factors, such as reaction rate, 
in the origmal Tanford and Pease equation (reference 6). 
Inasmuch as the present data are self-consistent, it is valid 
to use these data to predict the effect of initial mixture 
temperature on maximum flame speed. 


SUMMARY OF RESULTS 


An investigation of the effect of initial mixture temperature 
on the laminar flame speeds of methane-air, propane-air, 
and ethylene-air mixtures gave the following results: 

1. The flame speeds of methane, propane, and ethylene 
increased with initial mixture temperature at an increasing 
rate. The percentage change in flame speed with change in 
initial temperature for the three fuels followed the decreasing 
order, methane, propane, and ethylene. Empirical equations 
for maximum flame speed u (cm/sec; based on the outer 
edge of the shadow cast by the cone) for use in the tempera- 
ture ranges indicated were: 
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For methane, at an initial mixture temperature T, of 141° 
to 615° K, 
u=8+0.000160 7, 


for propane, for 7,—200? to 616° K, 
u=10+0.000342 Tj» 

and for ethylene, for 7,—200? to 617° K, 
u=10-+0.00259 Ty! 


2. Flame-speed values based on the outer edge of the 
shadow cast by the Bunsen cone (total-area method) were 
close to those based on both visible Bunsen cone measure- 
ments and tube measurements for methane and propane, but 
the variation in values obtained by the three methods was 
greater for ethylene. For each fuel, a factor was determined 
to allow conversion of values based on the inner edge of the 


shadow to values based on the outer edge of the shadow; _ 


the conversion factors were constant with respect to both 
temperature and composition, within the accuracy of the 
measurements based on the outer edge of the shadow. The 
flame speed of propane-air flames was independent of tube 


diameter from 10 to 22 millimeters and of values of stream- _ 


flow Reynolds number from 1500 to 2100. | 

3. Both the thermal-theory equations presented by 
Semenov and the square-root law of Tanford and Pease 
(diffusion theory) were used to predict relative flame speeds 
within approximately 20 percent for the temperatures and 


—— Le. 


gases studied. The comparative agreement between the  _ 


values predicted by the different theoretical equations and 
the experimental value varied for the three fuels. In general, 
for elevated temperatures, the thermel-theory values were 
high and the square-root-law values were low; whereas, for 
low temperatures, the values predicted by the two theories 
were substantially the same. However, even qualitative 
observations regarding the relative merits of the two theories 
are subject to the assumptions made and the flame tempera- 
tures and activation energies used. 

4. Very good straight-line correlations between maximum 
flame speed and either the hydrogen-atom concentration 


alone or the summation of the effective relative concentra- 


tions of hydrogen atoms, hydroxyl radicals, and oxygen 
atoms were found for the three fuels. Equally good correla- 
tions were obtained either from flame temperatures based on 
sodium D-line measurements or from adiabatic flame 
temperatures. The slopes of the correlation lines were 
steeper than those reported elsewhere in investigations of 
the effects of mixture composition and fuel type on flame 
speed. Further experimental and theoretical investigation 
may show why the slopes were different and give better 
insight into the mechanisms involved in fame propagation. 
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